Quantum emulation of a spin system with topologically protected ground states using 

superconducting quantum circuits 
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Using superconducting quantum circuit elements, we propose an approach to experimentally 
construct a Kitaev lattice, which is an anisotropic spin model on a honeycomb lattice with three types 
of nearest-neighbor interactions and having topologically protected ground states. We study two 
particular parameter regimes to demonstrate both vortex and bond-state excitations. Our proposal 
outlines an experimentally realizable artificial many-body system that exhibits exotic topological 
properties. 
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I. INTRODUCTION 

Interesting phenomena, such as the Aharonov-Bohm 
effect and Berry phases, can occur in physical systems 
with nontrivial topology in real or parameter space. 
Topological quantum systems are now attracting con- 
siderable interest because of their fundamental impor- 
tance in diverse areas ranging from quantum field the- 
ory to semiconductor physics^ with the most recent ex- 
ample being the exploration of topological insulators^ 3 - 
These topological physical systems may also have po- 
tential applications because they are robust against lo- 
cal perturbations. Specifically, a topologically protected 
quantum state degeneracy cannot be lifted by any lo- 
cal interactions fiii It is therefore natural to consider us- 
ing topological phases for applications requiring a high 
degree of quantum coherence^ For example, it has re- 
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ccntly been pointed out that non-Abelian anyon 
a fractional quantum Hall system can lead to topologi- 
cal quantum computing^ Anyons are neither bosons nor 
fermions, but obey anyonic braiding statistics.— Unfor- 
tunately, they have not yet been convincingly observed 
experimentally in any physical system. 

Instead of only looking for naturally existing topolog- 
ical phases, one could also design artificial lattice struc- 
tures that possess desired topological phases. One ex- 
ample is the Kitaev honeycomb model^ which requires 
that the spin (natural or artificial) at each node of a 
honeycomb lattice interacts with its three nearest neigh- 
bors through three different interactions: <J x a x , cr v <j y , 
and a z a z . Depending on the bond parameters, this 
anisotropic spin model supports both Abelian and non- 
Abelian anyons^ Its realization could potentially lead to 
experimental demonstration of anyons and implementa- 
tion of topological quantum computing. However, the 
requirement for anisotropic interactions is tremendously 
demanding and generally cannot be satisfied by natural 
spin lattices. 

Various artificial lattices may possess interesting topo- 
logical phases. For instance, it has been proposed 



that a triangular Joscphson junction array may have a 
two-fold degenerate ground state that is topologically 
protected. 1 ^ 11 A recent proposal suggests the use of ca- 
pacitively coupled Josephson junction arrays to simu- 
late a two-component fermion model that has topological 
excitations*^ There is also a suggestion that a Joseph- 
son junction array with properly designed interactions 
and topology can be local-noise resistant pi 3 - With respect 
to the physical realization of the Kitaev model, there are 
proposals using neutral atoms in optical lattices i 14 ' 15 ! 16 
One similarity among all of these proposals, whether 
based on Joscphson junction arrays or on optical lat- 
tices, is that they all require extremely low temperatures 
and precise single-atom manipulations. The reason is 
that topologically interesting properties are not gener- 
ally contained in the symmetry of the system Hamilto- 
nian. Instead they are only emergent properties at very 
low temperatures. 



Here we propose a quantum emulation of the Kitaev 
lattice using superconducting quantum circuits (Ref. Qj] 
gives a brief summary of this work). As for the topic 
of quantum analog simulations, see Ref. [TH for a brief 
review. In our superconducting network, a Joscphson 
charge qubit is placed at each node of a honeycomb lat- 
tice. These charge qubits behave like artificial spins and 
are tunable via external fields i 19 ' 20 i 21 Each charge qubit 
interacts with its three nearest neighbors through three 
different types of circuit elements. One advantage of our 
proposal is that some circuit elements involved and their 
functionalities at low energies have already been demon- 
strated experimentally — for example, the <j z <j z and a x a x 
couplings between charge qubits have been studied in 



experiments 
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Here we show theoretically that they 



can indeed provide the needed anisotropic interactions 
when included in a honeycomb lattice. We then identify 
the ground states of this network in two different param- 
eter regimes and show that it can have both vortex and 
bond-state excitations. We also describe how they can 
be generated using spin-pair operations. 
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FIG. 1: (Color online) (a) Schematic diagram of the basic 
building block of a Kitaev lattice, consisting of four supercon- 
ducting charge qubits (labelled 1 to 4) : (i) Qubits 1 and 2 are 
inductively coupled via a mutual inductance M; (ii) qubits 
1 and 3 are coupled via an LC oscillator; and (iii) qubits 
1 and 4 are capacitively coupled via a mutual capacitance 
Cm- Inset: The three types of inter-qubit couplings are de- 
noted as x-, and y- or z-type bonds. Here each charge qubit 
consists of a Cooper-pair box (green dot) that is linked to 
a superconducting ring via two identical Josephson junctions 
(each with coupling energy E.j and capacitance Cj), to form 
a SQUID loop. Also, each qubit is controlled by both a volt- 
age Vi (applied to the qubit via the gate capacitance C g ) and 
a magnetic flux (piercing the SQUID loop), (b) A partial 
Kitaev lattice (honeycomb lattice) constructed by repeating 
the building block in (a), where a charge qubit is placed at 
each site. A plaquette is defined as a hexagon in the lattice. 
The plaquette operator is defined as W p = o\o\o\o%o\o% 
and is shown for a given plaquette p. 



II. KITAEV LATTICE BASED ON 
SUPERCONDUCTING QUANTUM CIRCUITS 

At low energies, superconducting (SC) qubits can be- 
have as artificial spins. Among the varieties of SC qubits 
(charge, flux, phase ) 19 i 20 i 21 and other hybrid a 24 i 25 ' 26 ), 
only charge qubits are known to interact with each 
other in all the individual forms of a x a x , & y <J y , and 
a z a z (via a mutual inductance, an LC oscillator, and 
a capacitance, respectively) i 22 i 27 ' 28 i 29 Therefore, to emu- 
late a Kitaev honeycomb lattice, we propose to build a 



two-dimensional SC circuit network based on SC charge 
qubits. More specifically, on a honeycomb lattice a 
charge qubit is placed on each node [Fig. [ljb)] , and one 
of the three circuit elements is inserted along each bond 
of the lattice (denoted as the x- 1 y-, or z-type bond). 
A building block of this lattice is shown in Fig. [lja), 
which consists of four charge qubits that are connected 
via an x-, ay-, and a z-type bond. Each charge qubit is 
a Cooper-pair box connected to a superconducting ring 
by two identical Josephson junctions to give it tunabil- 
ity: Each qubit is controlled by both the magnetic flux 
§i piercing the SQUID loop and the voltage Vi applied 
via the gate capacitance C g . 

Naively, a circuit element should maintain its basic 
characteristics when inserted in a larger network, at least 
in the linear regime. However, as it has been shown in 
previous studies of hybrid qubits ; 24 i 25 i 26 a superconduct- 
ing qubit based on one particular variable (for example, 
charge) can acquire characters of another (for example, 
flux) when additional circuit elements are added to it. 
Therefore, here we first clarify whether the different cir- 
cuit elements in our honeycomb network maintain their 
basic individual characteristics (particularly the forms 
and strengths of the interactions) at low energies when 
lumped together. 

We first write down the Lagrangian of the quantum cir- 
cuits, choosing the average phase drop across the two 
Josephson junctions of each charge qubit as the canonical 
coordinates. After identifying the corresponding canon- 
ical momenta, we then derive (this derivation is shown 
in the appendix) the total Hamiltonian of the quantum 
circuits as 

i x — link y — link 

+ K ^,k). (1) 

z — link 

Here the free Hamiltonian of the ith charge qubit is 

Hi = E c {m - n gi ) 2 - Eji(®i) cos(p h (2) 

where E c = 2e 2 Cs/A is the charging energy of the 
Cooper pair box, with the total capacitance Cs = 
2Cj + C g + C m , and A = C| — C* ; n t = -id/dip, 
the number operator of the Cooper pairs in the zth 
box (which is conjugate to (pi); n g i — C g Vi/2e the re- 
duced offset charge induced by the gate voltage and 
Eji{$i) = 2Ej cos(7r$i/$ ) the effective Josephson cou- 
pling energy of the ith charge qubit, with $o = h/2e the 
flux quantum. 

The three nearest-neighbor couplings, shown as the x, 
y, and z bonds in Fig. 1(a), are given by 

K x (l,2) = Mhh, 

K y (l,3) = -4^ 1 ($ 1 )Sj 3 ($3)sin^ 1 sin^3, (3) 
K Z (1,A) = E m {m - %i)(n 4 - rigi), 
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where 



where 



E m = 
h = - 



nCz(C g + C m ) 
A$ 



4e 2 C m 
A ' 

I -j— | COS ( 



(4) 



Here I c = 2irEj/$ is the critical current through the 
Josephson junctions of the charge qubits (we assume 
identical junctions for simplicity), while Ii is the circu- 
lating supercurrent in the SQUID loop of the ith charge 
qubit. Note that the coupling strength between nodes 1 
and 3 (along a y-link), £ oc (C g + C m ) 2 , is affected by the 
mutual capacitance C m that connects qubit 1 (3) with its 
nearest-neighbor along a z-link. Compared to the case of 
two qubits coupled by an LC oscillator^ where £ oc Cg, 
the capacitive inter-node coupling along the z-link in the 
present circuit greatly increases the inter-node coupling 
along the y-link because usually C m 3> C g . This is an im- 
portant and positive consequence when multiple circuit 
elements are introduced to create different inter-node in- 
teractions. 

At low temperatures, only the lowest-energy states of a 
superconducting circuit element are involved in the sys- 
tem dynamics, which is quantum mechanical. For the 
particular case of a charge qubit, where E c ^> Ej, the 
lowest-energy eigenstates are mixtures of having zero and 
one Cooper pair in the box, when the gate voltage Vi is 



near the optimal point e/C g (i.e. 



'</* 



i). Defining 



|0)j and as the two charge states having zero and 
one extra Cooper pair in the box, we now have a two- 
level system as a quantum bit, or qubit. In the spin-i 
representation based on these charge states |0)j = |t)i 
and \l)i = (i is the index of the nodes), the system 
variables can be expressed as 



U 



COS CO; = -a- , 



sm ipi 



"2 ff < 



(5) 



(i.e., 



Here we consider the simple case with n 
all gate voltages on the different nodes are identical: 
Vi = V g ) and = $ e for all qubits. The low-energy 
Hamiltonian of the system is then reduced to 

H = J X J2 v^k+Jy E CT K + J - E ° z i< 

a: — link y— link z— link 



(0) 



J X 

Jy 

Jz 

h z 



7T(f> f 



MT C sin z y-^- 1 
Z[Ej($ e )} 2 < 0, 



1 



\E m > 0, 
E c + 
\Ej(<S> e ), 



(7) 



with Ej($ e ) = 2i?jcos(7r$ e /$ ). The reduced Hamilto- 
nian ^ is the Kitacv model with an effective magnetic 
field with z- and x-componcnts. Here h x and h z play the 
role of a "magnetic" field. Since J y oc h 2 x , to maintain 
finite inter-qubit couplings, h x cannot vanish. Therefore 
our Hamiltonian represents a Kitaev model in an always- 
finite magnetic field, although the field direction can be 
adjusted. This Hamiltonian has an extremely complex 
quantum phase diagram because of all the (experimen- 
tally) adjustable parameters. Here we arc particularly in- 
terested in whether it has topologically-intcrcsting phases 
and when such topological properties might emerge. 



III. VORTEX AND BOND-STATE 
EXCITATIONS 

Below we focus on two particular parameter regimes 
of the finite-field Kitaev model of (O, under the general 
condition that the z-bond interaction dominates over the 
other interactions. In particular, when J z 3> J x , | J y \ 2> 
\h z \, \ h x \, we identify a vortex state excitation. This case 
is described in Section A below. When h z = but h x 
is of the same order as J x and J y , we identify a new 
excitation that we call the bond state. We describe this 
case in Sec. III.B. The vortex state is a known topological 
excitation in the zero-field Kitaev model, while the bond 
state is specific to the finite-field Kitacv model. 



A. Kitaev lattice with dominant z-bonds in a weak 
"magnetic" field 



We first consider the case when 

Jz > Jx, \Jy\ > \h x \, \h x \, 



(8) 



and treat V = J2i(h z a* + h x crf) as the perturbation. 
Using perturbation theory in the Green function formal- 
ism^ one can construct an effective Hamiltonian for the 
lattice: 



2hl 



z— link x— link 



(9) 



where A.£ z t x ) is the excitation energy of the state 
°i |ffo)j i- c -> the energy difference between states 
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Hamiltonian of the Kitacv lattice takes the form 



FIG. 2: (Color online) (a) Two types of vortex excitations wi 
and W2- A pair of vortices are generated along the horizontal 
direction for Wi (vertical direction for W2) by the spin-pair 
operator <r z = a z I (a v = a y a x ) acting on a z-link. (b) Two 
bond-state excitations bi and 62, which are also generated by 
the spin-pair operators a z and a y on a z-link. 



l#o) and \go). Here \g$) is the ground state of 
the unperturbed Hamiltonian, i.e., Hamiltonian (JB]) with 
the perturbation term V excluded. Note that the ef- 
fective Hamiltonian H' only contains contributions from 
the second-order terms because both the first- and third- 
order terms vanish. Including the zeroth-order term (un- 
perturbed Hamiltonian), the total Hamiltonian of the 
system can be written as 



h = j' x e °i° x k + Jy E °H + Jl * E 

a; — link y — link z — link 

(10) 

where the effective z- and x-couplings are 



V = J - 

° z — ,J z 

J 1 - J - 



2% 
Ae z 



(11) 



Below we focus on the Abelian excitations. When 
J'z ^ J'xi \Jy\i the dominant part of the Hamiltonian H 
is that along the vertical links, 



Ha = J z E a j a k- 



(12) 



2 — link 



Under Hq, the two spins along each z-link tend to be 
aligned opposite to each other (| Tl) or I IT)) m order 
to lower their energies. Indeed, the highly degenerate 
ground state | g) of H is an arbitrary vector in the Hilbcrt 
subspace spanned by \<?<j)i, where N denotes the 

total number of z-links and a = f , J.. Within the ground- 
state subspace of Hq and up to fourth order^ the effective 



(13) 



where 



J' 2 J 2 

l&J'l 

C r r cr 2 cr 3 cr 4 (T 5 cr 6 



t2 t2 

u x J y 

16JF : 



(14) 



Here W p is the plaquette operator for a given plaquette 
p [see Fig. QJb)]. The operator W p for any plaquette p 
commutes with the unperturbed Hamiltonian H : 



[H ,W p } = 0; 



(15) 



so that [H , H c g] = as well, and the ground states \g} w 
of H e g form a subset of the degenerate ground states \g) 
of Hq. It is straightforward to show that W 2 \g) = \g), or 
Wp|g) = ±|5). Since J c ft > 0, to minimize the energy of 
\g) w , we need 



W p \g) w = \g) l 



(16) 



In other words, the eigenvalues of the W p operators in 
the ground state \g) w are w p = 1 for all plaquettes p. 

When some plaquettes undergo transformations that 
lead to w p = 1 — ► —1, the system gets into an excited 
state. The lowest-energy excitation corresponds to the 
generation of a pair of vortices when w p = 1 — > — 1 for 
two neighboring plaquettes. In this excitation process, 
each of the two neighboring plaquettes acquires a phase 
7r, which is equivalent to the addition of a flux quan- 
tum $0 through each plaquette. As shown in Fig. 2(a), 
such an excitation can be generated by applying either of 
the following two spin-pair operators on the ground state 

\g)w- 



\z<) 



with 



(17) 



(18) 



Here the two operators erf (of) and Ii (erf) act on the 
ground state 1^)^ at the bottom and top sites of the zth 
z-link, respectively. This pair of vortices, generated by 
either erf or of, are topological states with an excitation 
energy of 



As = 4 J eff 



(19) 



above the ground state. As shown in Refs. and [Irl 
these excitations exhibit the braiding statistics of Abelian 
anyons. The ratio between this excitation gap for the 
anyons and J z is 



Ae 



JxJy 



< 1. 



(20) 
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For example, if J z ~ 10GHz and both J x and J y are one 
tenth of J z , this gap would be about 1 MHz, correspond- 
ing to a temperature of 0.1 mK. This small gap requires 
an extremely low experimental temperature for suppress- 
ing the thermal activation of the ground state to the vor- 
tex states. Note that a different perturbative approaches 
shows that in the parameter region J' z > J' x , \ J y | , the 
spin-pair operators of and a v - generally create both vor- 
tex and fermionic excitations. However, in the limit 
of J' z 3> J' x ,\J y \, the dominant excitations are vortex 
states*^ which is consistent with the conclusion drawn 
above. 



B. Kitaev lattice with dominant z-bonds in a 
uniform "magnetic" field along the x-direction 

If we stay in the regime where the z-bond couplings are 
dominant (J z 3> J x ,\J y \), but place each charge qubit 
at the optimal point where n g = i, so that h z — 0, a 
different quantum phase arises when \h x \ is comparable 
to J x , | J y \. In other words, we now consider the regime 

J*> J x ,\J y \,\h x \, and/i z = 0. (21) 

Here the zeroth-order Hamiltonian is again the coupling 
along the z-bonds: H = J z J2 z -imk a j a k ( n °tice that 
here the coupling strength is J z , not J' z \ with the same 
highly degenerate ground state \g) as discussed in the 
previous subsection. To clarify the low-energy excita- 
tion spectrum in this regime, we again use perturbation 
theory in the Green's function formalism to remove the 
linear terms and derive an effective Hamiltonian in the 
ground state sub-Hilbert space of Hq. Up to second- 
order, the effective Hamiltonian takes the form 

tfiff = -^eff £ 0*0%, (22) 
2 — link 

where 

K cS = ^. (23) 

<J z 

The spin pair operator P z = crjcf at a z-bond (again the 
two Pauli operators act on the bottom and top nodes of 
the particular z-bond) commutes with the zeroth-order 
Hamiltonian [P Z ,H ] = (although it anti-commutcs 
with the four plaquette operators W p connected to this 
z-bond). Similar to W p , the pair operator P z also has 
two eigenvalues p z — ±1. Thus the ground state \g)b 

of -ffgff should satisfy p z = 1 for all the z-bonds in the 
system. In other words, 

P,\g)b = \g)b, (24) 

for all z-bonds. Since no two z-bonds share a node in the 
honeycomb lattice, and the lattice is completely covered 



by all the z-bonds, we can solve the eigenstates of P z of 

(z) 

each z-bond and obtain the ground state of H^ s as 
1 N 

lff)6 = 2Jv7a®(IU)i + liT><)- (25) 
»=i 

This is a nondegenerate ground state, which forms a sim- 
ple subset of the highly degenerate ground states \g) of 
Ho. It is maximally entangled within each z-bond, but 
not entangled at all between different z-bonds. In other 
words, the two-spin correlation function decays to iden- 
tically zero beyond a z-bond. The lattice is now an en- 
semble of maximally entangled "spin" pairs that are com- 
pletely independent from each other. This ground state 
is reminiscent of (and simpler than) the dimerized va- 
lence bond solid state discussed in the context of spin 
Hamiltonians i 31 i 32 There valence bond states refer to a 
singlet | tl — IT) f° r the electron spins, which is dictated 
by the Coulomb interaction and Pauli principle between 
electrons. 

When the pair operators of and of are separately ap- 
plied to the ground state at the ith z-bond [sec Fig. 2(b)], 
the excited states 

\Zi}=a!\g) b , \Y$=S» i \g) b (26) 

are called a bond state — while the pair operators are dif- 
ferent, the states they generate are only different by an 
overall phase because \g)b is a factored state for all z- 
bonds. A bond state at the ith z-bond corresponds to 
the change of p z = 1 — > — 1 at that particular bond. It 
is 2K c g above the ground state in energy. Notice that a 
bond state is an excitation that is completely localized to 
a particular z-bond. Furthermore, bond states are gen- 
erated by the same pair operators that generate the vor- 
tex excitations, although the ground states of the system 
are different in these two cases. In contrast to \g) w , the 
ground state \g) b is nondegenerate, and the bond state ex- 
citations are very different from the vortex states. This 
transition from vortex excitations to bond states occurs 
when we vary the parameters of the system (i.e., tuning 
n g to 1/2 and reduce <E> e from close to $o/2 so that h x in- 
creases to the same magnitude as J x and/or J y ), during 
which the topological property of the system changes. 

IV. THE BRAIDING OF EXCITATIONS 

A vortex looping around another vortex can produce 
cither a sign change or no sign change to the wave func- 
tion. The first case is denoted as an e-type vortex loop- 
ing around an m-type vortex, and the second case cor- 
responds to an e-type vortex looping around an e-type 
vortex (see, e.g., Ref. [H for a more detailed discussion). 
This indicates anyonic statistics between the e and m 
vortex states. Therefore, braiding, which refers to mov- 
ing one quasi-particle around another, is an important 
tool to determine the statistics of the quasi-particles (in 
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FIG. 3: (Color online) Schematic diagram of the procedures 
for braiding excitations, (a) The operations Uh and U v for cre- 
ating excitations, which are achieved by successively applying 
spin-pair operators at z-bonds along the horizontal (Ph) and 
vertical (P v ) paths. Here the paths P v and Ph intersect at 
a z-bond. (b) A combined operation U h ~ 1 U~ 1 UhU v for both, 
braiding the excitations created in (a), and fusing them to the 
vacuum, (c) The operations Uh and U v for creating excita- 
tions, which are also achieved by successively applying spin- 
pair operators at z-bonds along Ph and P v , but the paths P v 
and Ph do not intersect at a z-bond. 



the present case the vortices). Here we show an alterna- 
tive procedure for braiding an excitation with another, 
which can be applied to both vortex and bond states. 

Let us consider two particular evolutions for the sys- 
tem. The first evolution U v contains spin-pair opera- 
tions of = of of applied to the ground state |<jf) at three 
successive z-bonds along the vertical path P Vl as shown 
in Fig. [21a). Here \g) = \g) w for the vortex case and 
\g) = \g)b for the bond-state case. The second evolution 
Uh contains spin-pair operations <5f = of/j applied at 
four successive z-bonds along the horizontal path Ph as 
shown in Fig.^a). After these two operations in series, 
the state of the system is UhU v \g), where 



TT ~z~z~z~z 

U h = a i a 3 a 2 a 1 , 



Ur 



(27) 



Now we turn the evolutions backward by applying U v 1 
and Ur 1 to the system successively, so as to fusei^ the 



excitations to the vacuum (i.e., the ground state) [see 
Fig. [2Kb)] . The final state of the system is now 



l*/> = u^u- l u h u v \ 



(28) 



When the paths P v and Ph intersect at a z-bond, such 
as in the example given in Fig. [3Ka), where a\ and trf 
anti-commute, Uh and U v anti-commute as well: UhU v = 
— U v Uh- The final state thus becomes 



\*f) = -\g) 



(29) 



In other words, a phase flip e t7T resulted from the evolu- 
tions. For vortex excitations, this is equivalent to the case 
of an e-type vortex looping around an m-type vortex, as 
shown in Ref. [l6l . In contrast, when similar operations 
are applied but the paths P v and Ph do not intersect at 
a z-bond [see Fig. (3Jc), for example], UhU v = U v Uh, so 
that 



l*/) = l.9), 



(30) 



yielding no phase flip in the final state as compared to 
the initial state. For vortex excitations, this is equivalent 
to the case of an e-type vortex looping around another 
e-type vortex. 

The braiding of excitations, i.e. whether there is or 
there is no phase flip, can be revealed by means of 
Ramsey-type interference ! 16 ' 33 To achieve this, one can 
keep the same U v as above, but use 



where 



(5?)* = (trf)'Ji 



(31) 



(32) 



i.e., each of is replaced by half of the rotation. In the 
braiding case shown in Fig. 02a), 



jI (33) 
at the crossing point of paths Ph and P v . Thus, 

I*/} = U^U- x U h U v \g) = (v z 2 )-jW z 2 rH~g) 

= i(a z 2 )- 1 \g)=ia z 2 \g)=i\Z 2 ), (34) 

similar to the case with an e vortex braiding with a super- 
position state of an m vortex and the vacuum^ However, 
in the case without braiding [see Fig.^c)], 



\*f)=U h - 1 U- 1 U h U v \g) = \g). 



(35) 



Therefore, the braiding of excitations can be distin- 
guished by verifying if an excited state \Z 2 ) occurs at 
the crossing point of paths Ph and P v . 

While the vortex state described by Eq. ([13)) and the 
bond state described by Eq. (f2"2"| are very different ex- 
citations, they have similar braiding properties. In the 
braiding procedure shown above, the system is initially 
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in the vacuum (cither \g) w or \g)b)\ after the braiding op- 
erations in Eq. ([25]) , the system is fused to the vacuum 
again, but with a sign change to the ground-state wave 
function no matter which ground state the system starts 
with. In order to distinguish the difference between the 
vortex and bond-state excitations, one needs to focus on 
the intermediate steps of the braiding operations. Take 
U v in Eq. 1(27]) as an example. When it is applied to \g) w , 
the spin-pair operator 5\ in it first creates a pair of e vor- 
tices, and then the other spin-pair operations o| and Og 
successively move one vortex downward along the vertical 
path P v . The final state U v \g) w is also a pair of vortices, 
but the two vortices are separated by three z-bonds in 
the vertical direction [see Fig. [3]Ja)]. Importantly, this 
pair of vortices U v \g) w is degenerate with the pair of vor- 
tices cr\\g) w - However, in sharp contrast to the vortex 
case, when U v in Eq. (|27|) is applied to | (?) & , each of the 
spin-pair operations of , i = 1,2,3 creates a bond state 
and the final state U v \g)b is nondegenerate with the bond 
state a\\g)b- 

V. IMPLEMENTATION OF QUANTUM 
ROTATIONS 



—icrf, so the of operation on the ith qubit is given by 

of = e^Rf (tt), (38) 

while half of the rotation is 

(of)3= e W 4 i?f(7r/2). (39) 

The corresponding inverse rotations can be achieved by 
shifting the gate voltage to n g i < n g . 

A a x rotation of a particular charge qubit can be gen- 
erated by a similar approach. Specifically, when n g = i 
and $ e = 0, one has h z = 0, h x = —Ej, and J x = 0. 
Again the honeycomb lattice is separated into a se- 
ries of one-dimensional chains. Here we assume that 
Ej | J y | , J z , achievable in this charge-qubit system, 
which allows us to perform a single-qubit rotation driven 
by Ej. When n g = i, for a time r we switch off the 
flux in the SQUID loop of the ith qubit (the working 
point of this Kitaev lattice is usually at < < $o/2), 
producing a local x-type rotation on the ith qubit: 

R* (0) = exp[i(SEjT/h)af] = exp(i0of /2), (40) 



As indicated in previous sections, single-qubit rota- 
tions are needed to create vortex and bond-state exci- 
tations, and to perform braiding operations. Below we 
show that these quantum rotations of individual qubits 
in the honeycomb lattice can be achieved via electrical 
and magnetic controls. The key is to reduce the coupling 
between a specific qubit and its neighboring qubits to 
such a degree that its single-qubit dynamics dominates 
for a period of time. 

To generate a <j z rotation at a particular charge qubit, 
we consider the following approach by controlling both 
the magnetic flux through SQUID loops and the local 
electric field. Specifically, when the magnetic flux in the 
SQUID loop of each charge qubit is set to <E> e = &o/2, 
h x — and J y = 0, so that the honeycomb lattice is now 
decoupled into a series of one-dimensional chains. For a 
charge qubit E c ^ Ej($> e ), thus E c ^> J x . We further 
assume that E c E m , so that E c J z as well. One can 
now shift the gate voltage for a period of time r at the 
ith lattice point far away from the usual working point 
n g ~ | of the Kitaev lattice, so that the corresponding 
single-qubit energy 5h z ~ E c (instead of ~ Ej) is much 
larger than both J z and J x . Such a parameter regime 
should be reasonably easy to achieve for a charge qubit. 
This operation of shifting n g should yield a local z-type 
rotation on the ith qubit: 

Rl{9) = exp[-i(6h z T/h)(rf} = cxp(-i6»o.f /2), (36) 

where 

Sh z = h z (n gl ) - h z (n g ). (37) 
When 9 = 26h z T/fo = tt (where n g i > n g ), R*:(tt) = 



where 

6Ej = Ej-^Ej($ e ). (41) 
The of rotation on the ith qubit is 

of = e -"/ 2 i?f (tt), (42) 

where 25Ejt/?i = tt. Note that when the flux in the 
SQUID loop of the ith qubit is switched off to produce 
a local cc-rotation, the flux in the SQUID loop of the 
nearest-neighbor qubit that is connected to the ith qubit 
via an LC oscillator should be simultaneously shifted to 
a value around $0/2, so as to keep \J y \ between these 
two qubits much smaller than Ej. 

With both of and of rotations available for the ith 
qubit, the a\ rotation is given by 

a v = e -*«/* (T * a f. (43) 

Therefore, one can construct all the wanted operations 
of and of for generating both vortex and bond-state ex- 
citations by using the single-qubit rotations of and of. 

In order to obtain accurate z- and x-type single-qubit 
rotations, we assume that E c and Ej are much larger 
than the inter-qubit coupling. Actually this somewhat 
stringent condition can be loosened for realistic systems. 
As shown in Ref. H3, accurate effective single-qubit ro- 
tations can still be achieved using techniques from nu- 
clear magnetic resonance when the inter-qubit coupling 
is small compared to single-qubit parameters (instead of 
much smaller than E c and Ej). 
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VI. DISCUSSION AND CONCLUSION 

In this paper, our main objective is to construct an 
experimentally feasible proposal to emulate the Kitaev 
model on a network made of superconducting nanocir- 
cuits. To focus on the topological properties of the sys- 
tem, we choose the limit of identical qubits and identi- 
cal coupling strength. Furthermore, we fix the mutual 
inductances and the capacitances of the various circuit 
elements involved. There are basically two tunable pa- 
rameters: the gate voltage on the Cooper pair boxes (n g ) 
for each charge qubit, and the magnetic flux <!> e through 
the SQUID loops connected to the Cooper pair boxes. 
Within the regime where z-bonds dominate in interac- 
tion energy scale ( J z much larger than all other couplings, 
including J x , J y , h Xl and h z ), we have explored two lim- 
iting cases: one with weak effective magnetic fields {\h x \ 1 
\h z \ <C J x , \J y \), the other with the effective field only 
along the x-direction. We have identified some proper- 
ties of the relevant ground states, and the low-energy 
excitations, the vortex and bond states. However, much 
more study is needed to completely clarify the energy 
spectrum, the phase diagram, and the dynamics of this 
superconducting network. 

One observation we have made is that the vortex exci- 
tations and bond-state excitations can be generated us- 
ing the same spin-pair operations, starting from different 
ground states (\g) w and \g)b) that depend on the system 
parameters. We have also shown that while \g) w is highly 
entangled, \g)b is only entangled locally but not globally. 
This quantum phase transition requires more extensive 
studies to identify the critical point and related critical 
phenomena, such as how system entanglement changes 
near the transition point, and most importantly how its 
topological properties change. It would also be worth- 
while to investigate the system spectrum (from vortex 
excitation to bond state excitation) and dynamics dur- 
ing this transition, similar to our study of quantum phase 
transitions between Abelian and non-Abelian phases of 
the Kitaev models While such studies are generally nu- 
merically intensive, it would help reveal the exotic topo- 
logical properties of this many-body model. 

With the elementary building blocks given in Fig.QJa), 
one can construct Kitaev spin models on other types of 
lattices as well (see, e.g., Refs. IH and[13). In particular, 
it has been shown that in the absence of a magnetic field, 
the Kitaev model on a decorated honeycomb lattice 3 -^ 
can support gapped non-Abelian anyons. The quantum 
analog simulation of Kitaev models on different lattices 
using superconducting circuits should shed light on the 
novel properties of these topological systems. 

There are two important open issues in the study of 
building a superconducting qubit network to emulate 
a spin lattice. One is the role played by the decoher- 
ence of individual qubits, and the other is the measure- 
ment of correlated states on a qubit network. It is well 
known that charge qubits suffer from fast decoherence. 
However, it is not clear how decoherence would affect 



the topological excitations. Indeed, the faster decoher- 
ence of charge qubits may allow the system to reach 
its ground state faster. Furthermore, topological exci- 
tations are supposed to be robust against local fluctua- 
tions, so that decoherence in individual nodes may not 
easily destroy excitations such as the vortex state. Quan- 
tum measurement is another open issue in the study of 
collective states, whether ground states or low-energy ex- 
citations, of a qubit lattice. While single-qubit measure- 
ment of superconducting qubits can now be done with 
quite high fidelity ) 38 i 39 and two-qubit correlation mea- 
surements have been done^ measuring multi-qubit cor- 
relations requires further theoretical and experimental 
studies. We hope that our proposal acts as another incen- 
tive for researchers in the field of superconducting qubits 
to look for ways to perform measurements that can reveal 
quantum correlations. 

In conclusion, we have proposed an approach to em- 
ulate the Kitaev model on a honeycomb lattice using 
superconducting quantum circuits, and shown that the 
low-energy dynamics of the superconducting network 
should follow a finite-field Kitaev model Hamiltonian. 
We analytically study two particular limits for system 
parameters, explore their ground state characteristics, 
and identify their low-energy excitations as vortex states 
and bond states. We further show that both vortex- 
and bond-state excitations can be generated using the 
same spin-pair operations, starting from different ground 
states. Our proposal points to an experimentally realiz- 
able many-body system for the quantum emulation of the 
Kitaev honeycomb spin model. 
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APPENDIX A: DERIVATION OF THE 
HAMILTONIAN 

Below we derive the Hamiltonian of the honeycomb 
lattice constructed with superconducting quantum cir- 
cuits as described in Fig. [TJ For simplicity all charge 
qubits have the same parameters. Furthermore, the mu- 
tual inductances M, the LC oscillators, and the mutual 
capacitances C m for the x-, y-, and z-couplings are also 
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identical, respectively. Since the self inductance of the 
SQUID loop in each charge qubit is small, the voltage 
drop across this loop inductance can be ignored as com- 
pared with the voltage drops across the Josephson junc- 
tions in the loop. Also, we assume that the capacitance 
of the LC oscillator is much larger than the gate capaci- 
tance and the mutual capacitance, i.e., C>C 9 , C m . The 
total electrical energy of the qubit lattice can be written 
as (the 1-2 ad 1-3 couplings are magnetic and will be 
discussed later) 



T 



(Al) 



-link 



where the summation is over all the z-links. the term 
Ti4 contains the charging energies of the nodes on either 
end of a z-link in the building block, together with the 
coupling across the link. It is given by 



14 



i=l,4 v 7 



Pi 



$0 



CgVgi 



'<!>o 
1 2^ 



_c$ 2 



Li 



(A2) 



where 



ai 

(X4 



(C g 



C m )$L 
C m )$L< 



(A3) 



and C E = 2Cj + C g + C m . Here (® /2n)<Pi = V Jl is the 
average voltage drop across the two Josephson junctions 
of the ith charge qubit, and = Vl = Vv) is the 
voltage drop across the LC oscillator connected to qubit 
1 (4). 



The Langrangian of the qubit lattice is 



C = T-U, 



(A4) 



where U is the total potential energy of the system, in- 
cluding Josephson coupling energy and magnetic energy 
in all the inductors in the network. To derive the sys- 
tem Hamiltonian, we choose ipi and $l as the canonical 
coordinates. The corresponding canonical momenta are 
thus 



Pi 



PL 



dL 

dtpi 
dL 



(A5) 



More explicitly, 

PI = ^ - Cm (^) <P* 

§j) {a x + C g V gl ), 

PL = C$L - Cm (J^J (tf4 + <fs) 

+(C g + O m )(^y i p 1 + i p 2 ), (A6) 

where the subscript 5 denotes the qubit which is con- 
nected to qubit 2 via the mutual capacitance C m . In the 
limit of C S> C g , C m , pl ~ C<I> l ■ Thus one has 

CsXi + C m XA 



<Pi 



(0>o/27r) 2 A ' 
C m Xi + C%X4 

($o/2tt) 2 A ' 
Pl_ 
C ' 



(A7) 



with A = C| - C m , and 



( ^ ) 



(A8) 



The Hamiltonian of the honeycomb lattice is thus 
H = ^ (pilfi +P4P4 +Pl$l) - C 



z-link 



= T 14 + ^ 



(A9) 



z— link 



where 

r li = ^§^ + ^l- + ^^ + §;. (A10) 

2(I^) 2 A 2(^) 2 A (£) a A 2C 

We now perform two gauge transformations, so that p\ 
and P4 become: 



(All) 



- [ — I &4 = P4 • 



After these gauge transformations, tpi and (p4 become 



C E /2tt 



C s /2tt 



(A12) 
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and 7i4 can be expressed as 



with 



T 14 = \c^K 2 {px,V g x) + l -C^JC 2 {p Al V gi ) 



vl 



+C m )C^ 1 ,V gl )K(n, V gi ) + ^, (A13) 



where 



(If) a 1 / 2 



(A14) 



Based on these building blocks, instead of the z-links, the 
Hamiltonian of the qubit lattice can now be rewritten as 



BB 



U. 



(A15) 



BB 



Here the summation is over all the building blocks and 
7bb, for a building block shown in Fig. HJa), is given by 

T BB = ic s /C 2 (p!, V g i) + ic s /C 2 (p 2 , V g2 ) 

+ ^s/C 2 (p 3 , V g3 ) + ic E /c 2 (p 4 , V g4 ) 



-C m K.(j>l, Vgl)JC(p4, Vgi) 



2C" 



(A16) 



where the prefactor | in the second, third and fourth 
terms (instead of | as in the first term) is due to the 
lattice geometry that each of the qubits 2-4 is shared by 
three building blocks. 

The potential energy of the system consists of the 
Josephson energy —.Eji (<5il) cos fi of each qubit, the 
magnetic energy $|/2L of each LC oscillator, the self- 
inductance energy \L q If of each qubit, and the mutual- 
inductance energy —Mlilj between every pair of nearest- 
neighbor qubits coupled via M. In particular, the 
Josephson coupling energy is 



E^u.) = 2Ejcos 



(A17) 



The supercurrent in the SQUID loop of the ith qubit is 
= — I r sin ( " ~ ^ ) cosifi . (A18) 



$ 



Here I c = 2ttEj/^q is the critical current of the Joseph- 
son junction, $0 = h/2e is the flux quantum, L q is the 
SQUID loop inductance of each qubit, and the total mag- 
netic flux in the loop of qubit i is given by 



$ jL = Qi + Lqli-MIj, 



(A19) 



with $j the externally applied magnetic flux in the loop 
of qubit i and Ij is the supercurrent in the loop of qubit 
j that is coupled to qubit i via M. Based on the building 
blocks, the potential energy U can been written as 



U = U BB, 
BB 



(A20) 



Ubb 



-Ejx(&il)cOS(Pi - -Ej 2 (<&2l) COS f 2 

~E J3 ($ 3L )cOSf 3 - ^E J4 ($ 4L )cOSf 4 
$2 



1 



1 



2L + 2 Lqh + 6 Lqh 

~L q ll-Mhh, 
6 



-L n 



9^3 



(A21) 



where the prefactors ^ and ^ are again due to the lattice 
geometry that each of the qubits 2-4 is shared by three 
building blocks. 

Usually, the self inductance L q and the mutual induc- 
tance are much smaller than the Josephson inductance 
of each junction in the qubit loop. Thus, one can expand 
Eqs. (|A17|) and (|A18|) around ir<& l /$ and keep the lead- 
ing terms, as in Ref. The potential energy is then 
reduced to 

Ubb = -Ejx{$x) cos fx - ^E J2 ($ 2 ) cos^ 2 

--E J3 ($ 3 )costp 3 - \-E J4 (§ 4 )cosf A 

$2 

H — - 
2L 



Mhh, 



(A22) 



where the supercurrents are replaced by 



L = -I c sin 



$0 



cos tpi 



(A23) 



In Eq. (|A22|) . we have also omitted constant terms which 
are reduced to identity operators in the qubit subspace, 
because these terms only shift the zero energy of the sys- 
tem. 

Using the gauge transformation 



C E /2tt 



when the fluctuations of a,- are weak so that2S 



one has 

cos if 



COS (fi 



C s /2tt 



A V$ 

The potential energy Ubb is given by 
Ubb 



sin fi 



(A24) 



(A25) 



(A26) 



1 



■Ej\{$x) COS If 1 - -E J2 ($ 2 ) cos^ 2 
-3 Ej 3 (® 3 ) COS f 3 - ^E J4 (<P 4 )cosp 4 



1 

2L 



<I>/ 



/2ttL 



{Yx + Y 2 ) 



2ir 2 L 



{Y 1 +Y 2 f + MhI 2 , 



(A27) 
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with 



Yi 



Cs(Cg + C„ 

A 



(A28) 



Here the terms modifying the Joscphson coupling en- 
ergy are ignored because they are much smaller than the 
Josephson coupling energy. 

The term p 2 L /2C in Eq. (|AT6|) is the ki- 
netic energy of the LC oscillator and the term 



1 

2L 



<!»/ 



^r)(Y 1 + Y 2 



in Eq. (|A27|) is the poten- 
tial energy of the LC oscillator. When the frequency 
of the LC oscillator is much larger than the qubit 
frequency, the LC oscillator remains in the ground 
state, so that these terms can be removed from the 
Hamiltonian^i Thus, the Hamiltonian of the system can 
finally be written as 



i x— link y— link 

+ K ^,k). (A29) 



z — link 



Hero 



Hi = -^^{p^Vgi) - Eji{$i) cos&, 



(A30) 



with K (pi, V g i) given in Eq. (j A14|) . For the building block 
shown in Fig. QJa), the three nearest-neighbor couplings 
K x , K y and K z are given by 

K x (l,2) = Mhh, 

#2,(1,3) = -4^ J1 ($i)Sj 3 ($3)sin^i sin^ 3 , 
K y (l,A) = C m /C(p lj F ff i)/C(p 4 ,'M, (A31) 



where 



Z = L 



nC-zjCg + C m ) 
A$ 



(A32) 



In Eq. (|A29|) . the terms with sin 2 ^ are also removed 
because they are reduced to the identity operators in 
the qubit subspace. The canonical coordinates ifi and 
momenta pi are conjugate variables, and they obey the 
commutation relation: 



Below we give two examples of parameter regimes 
where the physics we discussed in this paper can be re- 
alized. For a quantum circuit with two charge qubits 
coupled by a mutual capacitance, the typical parameters 
are Cj w 500 aF, C m « 30 aF, C g « 0.5 aF, and Ej w 
15 GHz (see, e.g., Ref.[H). Here we choose C m = 200 aF 
so as to have a stronger capacitive coupling, Cj w 400 aF, 
and C g w 0.5 aF. These parameters give E c ~ 80 GHz 
and J z ?s 8 GHz. We also choose Ej = 20 GHz and 
apply a magnetic flux $ e in each qubit loop such that 
Ej($ e ) « 4 GHz. This gives \h x \ « 2 GHz. Because h z 
can be independently controlled by the gate voltage, it is 
easy to obtain \h z \ ss \h x \. Finally, wechooseM w 6.6 nH 
and the parameters of the LC oscillator are chosen as 
L w 3.8 fiH and C = 4C m = 800 aF. We then have 
J x w \J y \ Rj 4 GHz. The parameter regime given in 
Sec. III. A (i.e., J z ^> J x ,\J y \ ^> \h z \, \h x \) can thus be ap- 
proximately achieved. Also, J z is much smaller than the 
frequency of the LC oscillator cj = 1/ \JLC w 20 GHz, so 
that the lattice dynamics can be reasonably described by 
the Kitaev model in this regime. Note that E c m 80 GHz 
and Ej = 20 GHz, which are much larger than J z . Thus, 
the local quantum rotations erf and af for generating 
topological excitations can also be achieved. Though E c 
and Ej are much larger than or comparable to the fre- 
quency of the LC oscillator, the local quantum rotations 
arc implemented by changing the external fields applied 
locally on the qubits involved. It is expected that the to- 
tal Kitaev lattice will not be affected so much by these lo- 
cal operations because the topological properties should 
be robust against local fluctuations. 

For the parameter regime of Sec. III.B, we choose L w 
2.5 /iH, M sa 4.9 nH, and n g = h. The applied magnetic 
flux in each qubit loop is such that Ej($> e ) = 3 GHz. 
Other system parameters are chosen to be the same as 
those in the case above. Thus, we have J z f» 8 GHz, J x m 
\Jy\ ~ | /is | ~ 3 GHz, and \h z \ = 0. These parameters are 
much smaller than the frequency of the LC oscillator 
uo ~ 20 GHz, allowing us to consider only the ground 
state of the oscillator. Thus the Kitaev lattice can also be 
realized in this regime. Moreover, because E c « 80 GHz 
and Ej = 20 GHz, which are much larger than J z , the 
local quantum rotations af and af at the ith site can be 
implemented. 



[<Pj,Pk]=ifiSjk, (A33) 

where pj = —ihd/dtpj. Defining n.; = pi/h, one obtains 
Eq. (fT]) by replacing rij and 0i in Eq. (|A29[) with m and 



With the parameters considered here, the vortex exci- 
tation energy would be of the order of 0.1 GHz or larger, 
corresponding to an experimental temperature of 10 mK 
or higher, already accessible by currently available dilu- 
tion refrigerators. 
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